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Abstract 

We have been monitoring M31 with HRC and ACIS 
onboard Chandra regularly in the past two years. By com- 
bining eight Chandra ACIS-I observations taken between 



that the unresolved X-ray emission in the core region is 
much softer than most of the resolved X-ray sources in 
that region. Fifteen X-ray point sources are newly associ- 



ated with globular clusters (Di Stefano et al. 2002), which 



1999 to 2001, we have identified 204 X-ray sources 



f 



w ithin wnerl combined with the previous ROSAT results (Sup 



17' x 17' region of M31, with a detection limit P er ct aL 2001 D brin S s the total number of M31 globular 



the cer 

of ~ 1.6 x 10 35 ergs _1 . Of these 204 sources, 21 are iden- 
tified with globular clusters, 2 with supernova remnants 
(one of them is spatially resolved with Chandra), and 8 
with p lanetary nebula . By comparing individual 



clusters with detected X-ray emission to 48. Chandra and 
XMM-Newton also discovered severa l bright (Lx > 10 3 7 



ergs 1 ) transients in M31 and M32 (Garcia et al. 2000a 



images , 

about 1 50% of the sources are variable in time scales of 
months. We also found 14 transients. Combining all the 
available transients found in literatures, there are 25 tran- 
sients in M31 and M32 detected by Chandra and XMM- 
Newton; we present some of the long-term lightcurves by 
using the HRC, ACIS and XMM-Newton data. The spec- 
tral shape of 12 sources is shown to be variable, suggest- 
ing that they went through state changes. The luminosity 
function of all the point sources is consistent with previous 
observations (a broken power-law with a luminosity break 
at 1.7 x 10 37 ergs -1 ). However, when the X-ray sources in 
different regions are considered separately, different lumi- 
nosity functions are obtained. This indicates that the star 
formation history might be different in different regions. 

Key words: galaxies: individual (M31) - X-rays: galaxies 



1. Introduction 

M31 was observed by Chandra and XMM-Newton soon 
after these observatories were launched. In the first obser- 
vation (8.8 ks) of the core region by Chandra in 1999 Oc- 
tober, 121 point sources were identified within the central 
17' x 17' region and the nucleus was nicely resolved into 
five point sources which were not seen in previous missions 
( parcia ct al. 2000a ). Moreover, a bright transient was dis- 
covered ~ 26" from the nucleus. A relatively deeper XMM- 
Newton observation (34.8 ks) was made in 2000 June; 116 
sources were detected down to a limiting luminosity of 
6 x 10 35 ergs" 1 (0.3-12 kpc; |Shirey et al. 2001| ) and a 

865 



|Carcia et al. 2000b; ^Jsbronc ct al. 200l|~ Shirey et al 



pulsating supersoft transient with a periodicity of 
s was dis covered (psbrone et al. 2001). Moreover, both 
Chandra ( Garcia et al. 2001a ; Primini ct al. 2000] ) and 
XMM-Newton ( Shirey et al. 2001 ) observations confirmed 



200 1| ; [Kong et al. 200l| ; parcia et al. 2001b| ). The bright- 
est of these reached a peak luminosity of Lx ~ 3 x 10 38 



ergs~ 1 ( |Kong ct al. 200l|) . 

We report herein a brief summary of our Chandra 
monitoring program of M31 in the past two years; it will 
mainly focus on the properties (identifications, temporal 
and spectral variability, and luminosity function) of point 
sources in the central ~ 17' x 17' region of M31 as deduced 
from eight separate ~ 5ks ACIS-I observations spanning 
~ 1.5 years. In addition, we will also present early results 
of an X-ray resolved supernova remnant and an updated 
"X-ray movie" of M31 up to 2002 January. 

2. Observations 

M31 was observed with Chandra regularly as part of the 
AO-1 and AO-2 GTO program during 1999-2001 (This 
program continues as GO observations in AO-3). The pro- 
gram is originally designed to search for transients in M31. 
These observations consist of a series of HRC snapshots 
(~ 1 ks) which cover the whole galaxy (see Figure |TJ); 
once a transient is found, a follow-up ACIS (~ 5 ks) ob- 
servation then images the newly-discovered transient, oth- 
erwise, ACIS points to the nuclear region. In this paper, 
we mainly focus on the ACIS-I (10, II, 12 and 13) data 
centered in the central 16.9' x 16.9' region of M31; this 
consists of 8 separate observations from 1999 October to 
2001 June, with exposure times ranging from 4 to 8.8 ks. 
The actual region covered by the Chandra observations is 
slightly larger than 16.9' x 16.9' due to different roll an- 
gles. Sources near the outer edge of ACIS are therefore 
not observable in all eight exposures. In addition, we also 
make use of the HRC data and other ACIS pointings to 
construct the long-term lightcurves of sources. 

In order to create a deep image suitable for the detec- 
tion of faint sources, the eight observations were combined 
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1993|). The remaining 127 (= 204 - 77) Chandra sources 




Figure 1. The field- of -view of HRC-I (blue) and ACIS-I 
(red) pointings overlaid on an optical Digitized Sky Survey 
image of M31. The aim-points are marked with a cross. 



were not detected in the ROS 'AT HRI catalog, presumably 
because they are below the ROS AT detection limit or are 
variable. This Chandra catalog extends ~ 5x fainter than 
the ROS AT HRI catalog. We identify 22 Chandra sources 
with globular clusters; twelve of these globular clusters are 
identified as X-ray sources for the first time. We find two 
matches with supernova remnant (SNR). In particular, 
one of the them is spatially resolved by Chandra (see § 5) . 
Eight planetary nebula (PN) are found to be associated 
with our Chandra sources. However, the X-ray luminos- 
ity of our PN candidates are at least 3 order of magnitud e 
brighter than that in our Galaxy (e.g. Kastner et aL 2001 ); 
the X-ray colors are generally harder as well. Therefore, 
it is possible that these PN are either very unusual or 
something other than PN. We suggest that these may be 
objects similar to GX 1+4 in our Galaxy (i.e. symbiotic 
stars with a neutron star companion). We searched for 
matches between stellar nova as listed in the IAUC con- 
temporaneous to and within the field-of-view of our Chan- 
dra observations, and found no matches within 3". We also 
searched for matches with OB associations, as O and B 
stars may be moderate X-ray sources (Lx < 10 33 erg s _1 , 
Bcrghoefer et al. 1997). While this is well below our detec- 
tion limit, a group of O and/or B stars may reach our de- 
tection threshold, and star forming regions could conceiv- 
ably harbor massive X-ray binaries. However, we found no 
matches within our search radius of 3". Five possible fore- 
ground objects are also found in our field based on their 
X-ray colors and coincidence of stars. 



into a single stacked image with a total integration time 
of 39.7 ks. Figure|| shows the stacked "true color" image 
of M31, which is a composite of images from soft (0.3- 
1.0 keV), middle (1-2 keV) and hard (2-7 keV) bands. 
Soft sources appear red, moderately hard sources appear 
green, and the hardest sources appear blue. The image has 
been corrected for exposure and smoothed slightly with a 
Gaussian (a = 0.5") in order to improve the appearance of 
point sources; diffuse emission around the nucleus is also 
clearly shown. Also shown in Figure || is the "true color" 
image of the central 2' x 2' region of M31 with 1/8" pixel 
resolution, with the possible nuclear counterpart (M31*) 
marked (Garcia 2001). 

3. Source Detection and Identification 
Discrete sources in the stacked image were found with 



WAVDETECT ( [Freeman et al. 2002[ ), a wavelet detection 
algorithm implemented within CI AO. The central 2' x 2' 
region was treated separately by using the 1 / 8" image. A 
total of 204 sources were detected above 2.5a. 

We cross-correlated the source list with existing cata- 
logs of M31 objects. We find that 77 Chand ra sources have 
counterparts in the ROS AT HRI catalog (Primini et al 



4. Temporal and Spectral Variability 

The eight Chandra ACIS-I observations described herein 
span nearly 2 years from 1999-2001. This is substantially 
longer than previous surveys by ROS AT (2 observations 
separated by ~ 1 year). In order to study long-term X-ray 
variability, we computed a variability parameter following 
Primini et al. (1993). We found 99 X-ray variables, cor- 



responding to ~ 50% of the total. By comparing ROSAT 



observations to Einstein observations 10 years earlier, Pri- 



mini et al. (1993) found that ~ 42% of the X-ray sources 



within central 7.5' region were variable. By comparing two 



XMM-Newton observaitions separated by six months Os 



brone et al. (2001) found that > 15% of the sources in the 



central 30' were variable. 

We have also discovered 13 bright transients for which 
the source is found in at least one observation with lu- 
minosity of >5 x 10 36 erg s _1 and is not detected in at 
least one of the other observations. Note that the luminos- 
ity limit covers typical outburst luminosities of soft X-ray 
transients and Be/X-ray binaries in our Galaxy. One im- 
portant transient was missed by this analysis because it 
had a peak luminosity below our "bright transient" thresh- 
old during the eight ACIS-I observations considered here. 
This object is XMMU J004234.1+411808 (= CXOM31 



Figure 2. Top: Stacked "true color" Chandra ACIS-I image (39.7 ks) of the central ~ 17' x 17' region of MSI. This 
image was constructed from the soft (red; 0.3-1 keV), medium (green; 1-2 keV) and hard (blue; 2-7 keV) energy 
bands. The pixel size is 1.96" and the image has been smoothed with a Gaussian (a = 0.5",) function. Bottom: Stacked 
"true color" Chandra ACIS-I image of the central ~ 2' x 2' region of M31. The pixel size is 0.123". M31* candidate 
is marked. 
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Figure 3. Lightcurves of 5 bright transients and globular cluster Bo 86 as seen in the past two years with Chandra 
and XMM-Newton. Observations from HRC-I (blue), ACIS-I (red), ACIS-S (green) and XMM-MOS/PN (yellow) are 
included. 
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Figure 4- Spectral changes of three of the sources. The first two from the left change from lower luminosity at softer 
state to higher luminosity with harder state. The bright transient in the right shows a harder spectrum during its decay 
from the peak of the outburst. 



J004234.3+411809), which psbronc et al. (2001)| suggest 
is an X-ray nova. By adding all the published X-ray tran- 



it al. 2001; Bhirey 2001; Kong et al. 2001 



sients in M31 found by Chandra and XMM-Newton (Os- 

|Gt 



0011); Trudolyubov et al. 2002[ ), we conclude that 



In addition to transients, we also found several sources 
which have interesting long-term lightcurve. For example, 
CXOM31 J004218. 5+411758 is in the M31 globular clus- 
ter Bo 86 and shows a possible ~ 200 day modulation. 
This is reminiscent of the Galactic source 4U 1820-30 in 
the globular cluster NGC 6624, which has a 176-d long- 
term modulation. It is worth noting that X-ray variability 
on long time-scales (from days to years) has been found 
in many low-m ass and high-mass X-r ay binaries in our 
Galaxy (see e.g. Kong fc Charles 200C ). 

remained in outburst for more than one year and finally We searched for spectral variability in all of the sources 
turned off in 2001 June (Figure 3). In addition, three HRC- using a method analogous to that described in Primini 
I data points indicate that the source underwent flarings et al. (1993~)| , but replacing the flux with hardness ratio. 



bronc 
et al. 

there arc 25 "bright" transients in M31 (and M32) in the 
past two years. 

Figure 3 shows five representative lightcurves of tran- 
sients found in M31 and M32. The transient discovered 



with Chandra in the first observation (Garcia et al. 2000a) 



or state transitions during the early stage of the outburst. 
During the whole outburst, the energy spectrum is con- 
sistent with a power-law with slope of ~ 1.5 (see inset 
of Figure 3). The remaining three transients show typical 
lightcurves with a fast rise followed by an exponential de- 
cay like those seen in our Galaxy (Chen et al. 1997). The 
characteristic decay time is about 30-50 days. CXOM31 
J004309. 7+41 1901 is a recurrent transient; the source was 
detected with Einstein, ROSAT HRI and ROSAT PSPC 



at lum inosities ranging from 1 to 8 xlO 37 ergs 1 . [Buppoi 



et al. ( £001)| noted variability between PSPC exposures of 
a factor of 5. Except for the first Chandra observations, 
the source underwent a state transition to the low state at 



Only 12 sources show spectral variability, corresponding 
to 6% of the total population. This is of course a lower 
limit, as small changes in the spectra of weak sources are 
undetectable with the number of counts accumulated in 
our 40 ks of merged data. In order to further investigate 
the nature of the spectral variations, we fit simple spec- 
tra to two of the brighter spectral variables. The fits show 
that as the counting rate increases, the spectrum becomes 
harder (see Figure 4). This is reminiscent of atoll and Z 
sources in our Galaxy (see e.g. Hasinger & van der Klis 
1989| ). The luminosity of these two sources ranges from 



Lx < 10 3e ergs 1 ; it re-brightened to 3 x 10 37 ergs nn 
2001 November and shows decay in 2002 January observa- 



(0.4 — 1.0) x 10 3S erg s _1 , which is higher than the typi 
cal luminosity (< 10 37 erg s _1 ) of atoll sources. However, 
this luminosity is similar to that of the Z sources, which 
are believed to reach the Eddington limit (Psaltis et al 



tions. Beginning in the summer of 2001, we also monitored 1995). The luminosity and spectral changes appear con 



newly discovered X-ray transients with HST to search for 
UV counterparts. This program allows us to follow the 
evolution of six transients in M31 during their outbursts in 



sistcnt with that seen in Z sources as they move along the 
"normal branch" . These are probably the first extragalac- 



both X-rays and UV; it will be the first step in extending 



the highly successful RXTEj ASM and optical follow-up 
of Galactic X-ray transients to our nearest neighbor spiral 
galaxy. 



tic Z sources to be identified, except for LMCX-2 (Smak 
Kuulkers 2000] ) . Continued monitoring may confirm the 



nature of these sources by allowing us to trace out the full 
Z-shape of the spectral variations expected if this analogy 
is correct. Spectral variation is also seen in transients dur- 



Figure 5. Left: "True color" Chandra ACIS-I image of CX0M31 J004327.7+411829. This image was constructed from 
the soft (red: 0.5-1 keV), medium (green: 1-2 keV) and hard (blue: 2-7 keV) energy bands. The pixel size is 0.496" and 
the image has been smoothed with a 0.496" FWHM Gaussian function. Right: Hardness-ratio (1-2 keV over 0.5-1 keV) 
image of the SNR, with X-ray contours overlaid. The contour levels are arbitrarily chosen (from 0.5-7 keV smoothed 
image) to best follow the surface brightness variations. Note the "hot spots" in the north-western arc and southern 
region of the SNR. 



ing their evolution away from the outburst. Figure 4 shows 
the spectral evolution of the brightest transien t (~ 3 x 10 38 
ergs -1 ) ever found in M31 ( Kong et al. 2001 ); the source 
becomes harder during the decay of the outburst. While 
the source is in its decay stage, it is still relatively bright in 
a recent observations taken in 2002 January (Lx <] 10 37 
ergs -1 ). 

5. Discovery of an X-ray Resolved SNR 

We have identified two SNRs in the central region of M31 
with Chandra, and both of them were detected by ROSAT 
previously. One of them (CXOM31 J004327.7+411829) is 
found to be extended in X-rays as a distinct ring-shaped 
object with a diameter of ~ 11" (33 pc). Figure 5 shows 
the "true color" X-ray image of the SNR; also shown is 
the hardness ratio image of the SNR. It is clear that bright 
spots can be seen in north-western and southern regions 
of the SNR. This SNR is also detected in several optical 



surveys (e.g. d'Odorico et al. 1980; Blair et al. 1981); it was 
identified as an irregular, faint and high [SII]/H Q (~ 1) 
SNR. The optical extension is about 10" — 18", with a 
crescent open in the south-eastern part of the SNR. Our 
X-ray image also shows similar morphology. 

Initial examination of a 40ks ACIS observations on 
2001 October indicates that a non-equilibrium ionization 
(NEI) model with interstellar absorption is appropriate 
for modeling the SNR's spectrum. It is also evident that 
the spectrum is dominated by a broad emission of O VIII 
line at 0.654 keV and a blend of Fe L shell lines and Ne 



K shell lines around 0.9 keV. The best fitting NEI model 
indicates a 0.3-7 keV emitted luminosity of ~ 5 x 10 36 
ergs -1 . 

6. Luminosity Function 

In Figure 6 we plot the cumulative luminosity function 
(CLF) for all detected sources in the stacked image, and 
also separately the LFs of inner bulge (region 1; central 
2' x 2'), outer bulge (region 2; central 8' x 8' excluding 
region 1), disk (region 3; central 17' x 17' excluding re- 
gions 1+2), and bulge (regions 1+2 combined). The LF 
for all sources has a break at ~ 2 x 10 37 ergs -1 , with 
otx = 0.38 before the break and a2 — 1.63 after the break. 
This result is in good agreement with previous ROSAT 



flPrimim et al. 1993[ ) and XMM-Newton flShircy et al 
[2001 ) measurements of the LF 



The CLF for the inner bulge is significantly different, 
with a break at a lower luminosity (~ 1.6 x 10 36 ) and a 
significantly flatter distribution at the faint end. We per- 
formed a two-sample Kolmogorov-Smirnov (K-S) test for 
the LFs of inner bulge and disk, and find that there is 
only 3% probability that they are drawn from the same 
distribution. Simulations show that the flattening of the 
CLF of inner bulge is not due to incompleteness. 

However, the counts in a CLF are not independent and 
therefore direct fits to it will underestimate the errors and 
may produce a biased best estimate of the slope. In or- 
der to more accurately estimate the errors and slopes, we 
used a maximum likelihood method (e.g. Crawford et al 
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Figure 6. Left: Luminosity functions for all sources (Re- 
gions 1+2+3) and bulge (Regions 1+2). Right: Luminosity 



for impulsive star formation. Within M31, we find that the 
LF of three regions we studied have breaks at a different 
luminosities. The inner bulge has this break at the low- 
est luminosity, and the luminosity of the break increases 
monotonically as we go out from the inner bulge. If the 
breaks do indicate the epochs of star formation events, 
then these events occurred most recently in the disk of 
M31 and further back in time as we move towards the 
nucleus of M31. 

As well as the monotonic shift in the break luminos- 
ity, there is a monotonic shift in the slopes of the LFs. As 
we move in towards the nucleus these slopes become pro- 
gressively flatter. This is somewhat difficult to understand 
in the context of the discussion above, because the most 
luminous sources would be expected to have the shortest 
lifetimes. Loss of these sources as they age would tend to 
steepen the luminosity function, but we find flatter lumi- 
nosity functions in the apparently older populations. 

It is interesting to compare the LF of these three re- 
gions of M31 to those of other galaxies. M31 is not the first 
galaxy found to show a break in its LF nor is it the first 



functi dns for inner bulge (Region 1), outer bulge (Region 
2) and disk (Region 3). 

1970) to determine the slopes in the differential luminosity 
functions (DLFs). For the inner bulge, the LF is roughly 
consistent with a single power-law with a ~ 0.7 although 
the appearance of the LFs (Figure 6) and the results of 
simulations indicate the presence of a brake to a flatter LF 
at the lowest fluxes. However, we note that the statistics 
in the inner bulge region are insufficient to confirm the 
presence and constrain the size of the break. The slope of 
the outer bulge and disk below the break is roughly con- 
sistent with 0.5-0.6. Further away from the inner bulge, a 
break appears near 10 37 ergs -1 and it shifts to higher lu- 
minosity when considering the disk region (see Figure 6). 
In addition, the slope becomes steeper (from ~ 1.0 to 2.2) 
when moving from the bulge to disk. 

The black holes and neutron stars that power many of 
the M31 X-ray sources have formed through the evolution 
of initially massive stars. Because of this, the X-ray LF 
traces the history of star formation and evolution of these 
massive stars in binary systems. Breaks in the LF may in- 
dicate indicate an impulsive star formation event. As the 
X-ray binaries age, their average luminosity shifts to lower 
values and therefore the location of the break may be an 
indication of the how long ago the star formation event 
occurred (Wu et al. 2002). Luminosity functions which 



to show different LFs in different regions; both M81 ( Ten 
nant et al. 2001| ) and M83 flSoria fc Wu 2002[ ) show similar 



do not show a break may indicate that star formation is 
still occurring. Chandra observations have measured the 
breaks in the LFs of several near by galaxies, e.g. M8 1 
( [Tcnnant ct al. 200 1| ), NGC 1553 ( [Blaiiton et al. 200l| ), 



behavior. In cases where a single slope is fit to the LF (ie, 
there is no clear break) the LFs of early type galaxies and 
the bulges of spirals tend to be steeper (a ~ 1.7) than 
those of spiral disks and galaxies with active star forming 
regions (a ~ 0.8, e.g, Prestwich 2001 ). The opposite seems 
to be the case within M31: the disk has a steeper LF than 
the bulge region. We speculate that this difference may 
be related to the location of the breaks in the M31 LF, 
which are at a lower luminosity that those seen in other 
galaxies. At these lower luminosities we may be sampling 
a different class of source, and the steepness of the LF may 
be due inclusion of this new class of faint sources rather 
than a loss of bright sources. We note that there is some 
evidence that we are sampling a different class of sources 
as we move out from the bulge, in that the fraction of 
sources which show variability decreases monotonically. If 
these sources have an intrinsically steeper LF than bright 
accreting binaries they may be able to cause the steep- 
ening of the LFs as we move from the bulge towards the 
disk. 



7. Conclusions 

In the past two years, Chandra and XMM-Newton have 
produced fruitful results on M31 (there are 7 papers re- 
lated to M31 in this meeting). Clearly, it is just a begin- 
ning while data are still being analyzed and will be coming 
in in near future. A series of papers about the point source 
properties, the diffuse emission, the X-ray emission from 
the central supermassive black hole, and supersoft source 
populations are being prepared by our Chandra team. On 



8 



the other hand, an extensive XMM-Newton survey on M31 
is undergoing and the results would allow us for a detailed 
studies of point source properties and populations on the 
whole galaxy. 
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